Localization and its consequences 
for quantum walk algorithms and quantum communication 
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The exponential speed-up of quantum walks on certain graphs, relative to classical particles 
diffusing on the same graph, is a striking observation. It has suggested the possibility of new 
fast quantum algorithms. We point out here that quantum mechanics can also lead, through the 
phenomenon of localization, to exponential suppression of motion on these graphs (even in the 
absence of decoherence) . In fact, for physical embodiments of graphs, this will be the generic 
behaviour. It also has implications for proposals for using spin networks, including spin chains, as 
quantum communication channels. 



The starting point for the field of quantum compu- 
tation is the realization that quantum algorithms may 
perform tasks more efficiently than known classical algo- 
rithms for the same problem The first example of 
a complexity separation was discovered by Deutsch and 
Jozsa [2j. Perhaps the most celebrated example is Shor's 
factoring algorithm 3], which is exponentially more effi- 
cient than any known classical algorithm. 

A third class of quantum algorithms stem from random 
walks, which have proved a fruitful tool for finding clas- 
sical algorithms . And this has led to investigations of 
whether quantum random walks might offer additional 
benefit for producing fast algorithms. The first realiza- 
tion was that quantum random walks in one dimension 
propagate at a rate that is linear in time, compared to 
the square-root rate for classical random walks |§] . This 
already suggests that quantum walks may offer possi- 
bilities for at least polynomial speed-up of algorithms. 
Thus the subsequent discovery of quantum walks that 
are exponentially faster than classical walks is particu- 
larly striking and suggestive 0, 0, E3] ■ 

The clear message is that quantum effects can enhance 
computation. In this paper we point out that the oppo- 
site is also possible and that quantum effects can cause 
exponential suppression of computations. 

Our focus is on certain continuous time walks (in con- 
trast to the discrete time "coined" walks mostly used in 
algorithmic applications) which, when carried out on ide- 
ally perfect quantum graphs, have exponentially faster 
propagation than a classical random walk on the same 
graph. We will apply well-established ideas from the the- 
ory of localization [9( to these quantum walks to show 
that when the graphs have imperfections, as they surely 
will in any real physical implementation, the propaga- 
tion of quantum information is suppressed exponentially 
in the amount of imperfection. 

Our conclusion is that quantum walks on physical 
graphs are not likely to be useful for algorithmic pur- 
poses. It is worth stressing at this point, however, that 
applications of quantum walks in quantum algorithms 
are thought of as being simulated walks in the memory 
of a quantum computer - as is indeed always the case 



in the classical applications where the underlying graphs 
are typically exponentially large and hence would require 
exponential resources to realize physically. 




FIG. 1: The Graph g 4 

We describe in detail a particular example of propaga- 
tion on a graph which, in the absence of imperfections, 
has the property that quantum evolution is exponentially 
faster than classical evolution; the effect of localization is 
particularly stark here. Later we point out the implica- 
tions for quantum information processing tasks on other 
networks. 

The particular case we will consider is the graph Q n , 
first studied in It consists of two joined branching 
trees with In + 1 vertical columns of nodes; Q4 is il- 
lustrated in Figure Q The question of interest is how 
rapidly a particle starting at the left-most node reaches 
the right-most node. Let us first describe the idea be- 
hind the exponential separation between classical diffu- 
sion and quantum motion on this graph (when the graph 
is perfect). Our treatment follows that in 0,0] closely. 

The model of classical diffusion is that the particle is 
equally likely to move from the node where it finds itself 
to any of the nodes to which it is linked. Explicitly we 
may take the motion to arise from the Markov transition 
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matrix 

! — 7 * 7^ j if nodes i and j connected 
i ^ j if nodes i and j not connected (1) 
da i = j 

where 7 is the constant jumping rate between two con- 
necting vertices, and di is the number of edges incident 
to the vertex i. The probability of being at node i at 
time t is governed by 

^ = -E^W(*)- (2) 

3 

A classical particle finding itself in the middle of the 
graph (not at the left- or right-most node or one of the 
nodes precisely in the centre) is twice as likely to move 
towards the centre as away from it. Thus a classical par- 
ticle starting at the left-most node will diffuse rapidly to 
the centre of the graph but then it will diffuse exponen- 
tially slowly from the centre of the graph to the rightmost 
node. 

We now analyze the quantum walk on Q n |6j start- 
ing in the state corresponding to the left-most node and 
evolving with the Hamiltonian given by 

(o|fl-|6) = M aM (3) 

where \b) represents the state of a particle at node b. 

The key observation is that although there are ex- 
ponentially many nodes, and a Hilbert space of expo- 
nential dimension, in fact the system only evolves in a 
(2n + l)-dimensional subspace. This subspace is spanned 
by states \j) (where < j < 2n), the uniform superposi- 
tion over all vertices in column j, that is, 

|^_ 2 -ain[i,2n-J]/a ^ (4) 
a£ column j 

In this basis, the non-zero matrix elements of H are 
(j\H\j±l) = -V2 7 

which is depicted in Figure as a quantum random walk 
on a line with 2n + 1 vertices. 

-s/2 -s/2 -V2 -V2 -V2 -y/2 -s/2 -V2 

• • • • • • • • • 
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FIG. 2: The elements of the Hamiltonian ff/7 for Q4 when reduced 
to a walk on a line |fj- 

In particular it is seen that (away from the left-most 
or right-most vertices) the quantum particle has equal 
amplitude to move to the left as to the right, and in 
particular it is not exponentially suppressed from moving 



right once it has reached the centre. (Observe that the 
central node, coming from the centre of the graph Q n , 
has a special status, but that will us not concern here.) 

We now turn to the main point of this paper. If we 
were to try and build this graph physically and get a 
quantum particle to evolve on it, it would be virtually 
impossible to do so perfectly. In particular, for exam- 
ple, the distances between the nodes will inevitably vary 
slightly from edge to edge. Our main observation will be 
that the theory of Anderson localization implies that 
this variability will lead to suppression of the quantum 
evolution so that in fact a quantum particle starting at 
the left of the graph can only travel a distance to the 
right proportional to an inverse power of the degree of 
variability. The probability of arriving at a point be- 
yond this "localization length" is exponentially small in 
the distance from the starting point. Thus if the right- 
most node is further away than the localization length, 
the particle is effectively prevented from reaching it. We 
emphasize that localization is very much a quantum (or 
more strictly a wave) phenomenon; it does not affect clas- 
sical particles diffusing. Indeed, Anderson localization is 
often characterized as the quantum suppression of clas- 
sical diffusion. We feel it is striking that this is a case 
where quantum effects are suppressing rather than en- 
hancing possibilities for information processing. 

The imperfections in the quantum evolution we con- 
sider are within the graph itself; namely the Hamiltonian 
varies a little from node to node. But the evolution is uni- 
tary. In particular, we are not concerned with any inter- 
action of the quantum system with the environment, in 
other words we assume that there is no decoherence . 
In fact we will consider a rather weak model of imperfec- 
tion for our graph. In the original (2™ +1 + 2" — 2)-node 
graph, one expects that the interaction connecting each 
pair of nodes will vary slightly from edge to edge. This 
breaking of the symmetry of the original system would 
have a substantial effect on the system, namely that the 
evolution would now no longer proceed simply through 
the states but leak into the entire Hilbert space. Our 
model of the variability of the Hamiltonian is less severe 
than this: we assume that the evolution still proceeds 
through the states but that this effective walk on 
the line is subject to a Hamiltonian that varies from node 
to node on the line. It is to be expected that the more 
general situation in which variability is allowed for all 
nodes/edges can only further suppress the quantum evo- 
lution. Thus our system is still considered to walk on a 
one-dimensional graph, but now our Hamiltonian is 

if e =ff + ]Te,|j)(j|, (6) 

j 

where H is the Hamiltonian J^J. The variability is in- 
troduced via the parameters ej. In the first instance the 
parameters ej will be taken independently from a Cauchy 
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FIG. 3: (Colour online) Propagation within the left side of Siooo 
at times i according to; (top) the perfect Hamiltonian eq. J3J> an d 
the disordered Hamiltonian eq. JfjJ with (middle) disorder 8 = 0.03, 
and (bottom) disorder <5 = 0.06. 



distribution with parameter 5: the density is 



P(e) 



1 



(7) 



A particular reason for doing this is that the Cauchy 
distribution is most amenable to analytic treatment. 
Later we will point out that taking the tj to be chosen 
from other distributions makes no qualitative difference 
to the conclusions (although it does lead to some inter- 
esting quantitative differences), ft will also be clear that 
we could encode variability in the Hamiltonian in various 
other ways than for example the off-diagonal terms 
could have random com pon ents. Many studies in the lo- 
calization literature [13, LL3j show that these possibilities 
do not make a qualitative change to our conclusion. 

All our conclusions are encapsulated in Figure [3J The 
first plot shows the flow of intensity with time for motion 
along left half of the columns of Q n (0 < j < n) accord- 
ing to the Hamiltonian (JSJ - the perfect case. The signal 
clearly travels smoothly (and indeed at constant speed) 
along the chain. The remaining two plots show what hap- 
pens when the Hamiltonian is now 10 with the tj taken 
from the Cauchy distribution with increasing amounts of 
variability, <5. In all cases, we numerically computed the 
probability |(j|V>(i))| of a particle being in the state 
|j) for seven times t, where \ip (t)) = exp(—itH/h) |5), 
7 = 1, H = 1, and n = 1000. 

When 5 is non-zero, the quantum wave packet travels 
a certain distance along the graph, but is then seen to 
stop. The distance that the packet gets to is reduced as 5 
increases (in fact, as we shall see the distance is inversely 
proportional to 5). 

This behaviour is a classic example of Anderson lo- 
calization 0- It arises since the eigenstates, \E a ), of 
the Hamiltonian © are exponentially localized in space, 
rather than being extended into space as in the case 



©. The modulus of the amplitude of each eigenstate 
is bounded by an exponential function 



\(E a \j) \ < N a e 



-j /lc 



(8) 



where l a is the localization length and N a is a normal- 
ization constant. It is simple to see why this behaviour 
of the eigenstates causes the time evolution witnessed in 
Figure The state of the system at time t may be writ- 
ten as 



h&(*))=X>a 



-itE a /h 



E a ) 



(9) 



where ip a = (E a \ip(0)) are the amplitudes of the initial 
state and E a are the energy levels. Now consider a quan- 
tum state initially concentrated at some point. This ini- 
tial state will only have substantial amplitude in energy 
eigenstates which are localized in regions close to that 
initial point. As the state evolves with time, the phase 
of these amplitudes can change, but not the modulus of 
the amplitudes. Thus for all time the state is expanded 
in terms of states which have small amplitude far from 
the initial point. 

For the particular model of Hamiltonian variability 0, 
we can explicitly describe the behaviour of the localiza- 
tion length of eigenstates as 6 varies, treating as negligible 
the influence of the anomalous end and centre points of 
the graph. Namely, it has then been shown 0] that 
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(10) 

where E a = E a — 37. It is not difficult to check that 
the largest localization length occurs for E a = 37, and 
for small S (i.e. S <C V87) it becomes Z max ~ y/Ej/S. 
Since the Hamiltonian has a random component, it will 
almost surely not have an eigenvalue equal to 37; but 
the localization length for E a = 37 is an upper bound 
for the localization length of any eigenstate. And hence 
this gives an estimate of the furthest distance a signal 
starting at a given point will get to. The plots in Figure 
13 bear this out. 

This then is the main conclusion of our investigation. 
A particularly quantum effect, Anderson localization, 
causes a particle undergoing a quantum walk in the graph 
to be effectively stopped from propagating from left to 
right. More precisely, consider a fixed value of the degree 
of variability 5, and a series of graphs Q n with n increas- 
ing. Once the length of the graph n is greater than about 
the maximum localization length ^/8 / f/S, the probabil- 
ity of the particle reaching the right-most node, having 
started at the left-most node is exponentially small in n. 

This is reminiscent of the behaviour of the classical ran- 
dom walk on the same graph. However it is fundamental 
that the exponents arise from quite different sources. In 
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particular we note that in the original problem of motion 
on the graph Q n , a classical particle diffusing from left 
to right travels rapidly from the left most node to the 
centre, and the exponential time it takes to get from left 
to right arises from the fact that it is difficult to get from 
the centre of the graph to the right-most node. In the 
case of a quantum particle on an imperfect graph, local- 
ization means that it can be exponentially difficult even 
to get from the left-most node to the centre of the graph. 

It is important to point out that localization is not a 
consequence of our choice of Cauchy distributed disor- 
der. It is a generic property of disordered systems, in 
one dimension, occurring for any distribution. For other 
choices (e.g Gaussian or uniform on an interval), the qual- 
itative features are the same as for the Cauchy case, but 
the quantitative details may differ; in particular, if P(e) 
has a second moment, the localization length typically 
scales like 5~ 2 as 5 — ► 0, rather than d^ 1 |12|. 

The graph Q n is of particular interest because of its 
relationship to quantum algorithms and the exponential 
separation between between classical and quantum prop- 
agation. However our results also have implications for 
other quantum information processing tasks. For exam- 
ple, there has been considerable interest recently in using 
spin chains and other networks for propagation quantum 
information . The message of this paper is that local- 
ization will present a fundamental difficulty for these pro- 
posals even in the absence of external interactions with 
the system. (See also [l6j|.1 

To end this discussion, we emphasise two differences 
of the model considered here to noisy quantum computa- 
tion. First, noise in quantum computers is usually mod- 
elled as sufficiently independent stochastic variations in 
the dynamics, in particular on a single memory qubit 
they vary with time; here, we have randomness in the 
description of the system itself, i.e. the Hamiltonian, but 
on the other hand the localization effect remains even if 
the particular random Hamiltonian is known. Second, 
we consider here free evolution of the quantum state on 
a spatially extended system; this means that any com- 
parison with quantum computers should not be with the 
circuit model (for which also techniques for fault-tolerant 
computation exist [17J), but rather with computation in a 
closed system, like quantum cellular automata |18j. It is 
not known whether fault-tolerant techniques for univer- 
sal quantum cellular automata in one dimension apply to 
the type of error considered in localization: fault-tolerant 
computing usually has to assume a degree of indepen- 
dence of errors, both spatially and in time, and although 
able to cope with certain correlations, the persistent "fail- 
ure" of an interaction in always the same way throughout 
the evolution seems to present a new challenge. 

Also, it should be pointed out that localization is 
known to occur for a single spin excitation moving in 
spin lattices, and it may be expected for small numbers 
of excitations; however, universal quantum computation 



happens in the regime of many excitations, where it is not 
known if localization presents an obstacle for information 
propagation. 

We remark finally that the systems we have considered 
here are different to stroboscopic (<5-kicked) models (e.g. 
those in [l^) which are closely related to the kicked rotor, 
where localization occurs for quantum chaotic reasons 
|2flj , and its consequences have been extensively ex plor ed 
by Shepelyansky and co-workers (see for example 21]). 
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